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1.  Introduction 


This  paper  develops  a  new  approach  for  the  study  of  generalized  stability  of  parametrized 
families  of  matrices  or  polynomials.  The  matrix  case  is  emphasized  in  the  exposition. 

Generalized  stability  [V]  of  a  matrix  (resp.  polynomial)  entails  that  its  eigenvalues 
(resp.  zeros)  lie  in  a  prespecified  domain  of  the  complex  plane.  The  classical  stability  re¬ 
quirements  result  upon  defining  the  domain  of  interest  as  the  open  left  half  complex  plane 
(continuous-time  case)  or  the  open  unit  disk  (discrete-time  case).  Practical  considerations 
relating  to  damping  ratio,  bandwidth,  vehicle  handling  qualities,  etc.,  are  commonly  ex¬ 
pressed  in  terms  of  the  generalized  stability  formulation,  with  respect  to  a  suitable  domain 
in  the  complex  plane. 

The  problem  of  obtaining  necessary  and  sufficient  conditions  for  the  Hurwitz  stability 
of  polytopes  of  matrices  or  polynomials  has  recently  been  considered  by  several  authors 
(see  for  instance  [BD]).  This  research  direction  was  inspired  largely  by  Kharitonov  [K], 
who  studied  Hurwitz  stability  of  families  of  real  polynomials  with  uncertain  coefficients. 
Specifically,  he  studied  “interval  polynomial”  families,  i.e.,  families  of  polynomials  with 
coefficients  belonging  to  fixed  compact  intervals  of  the  real  line.  He  showed  that  the 
Hurwitz  stability  of  such  a  family  of  polynomials  is  equivalent  to  that  of  four  distinguished 
“corner  polynomials.”  Bartlett,  Hollot  and  Lin  [BHL]  showed  that,  in  the  case  of  an 
arbitrary  polytope  of  polynomials,  it  suffices  to  check  the  edges  (The  Edge  Theorem). 
These  strong  results  do  not  hold  however  for  polytopes  of  matrices.  For  the  case  of  Hurwitz 
stability  of  the  convex  hull  of  two  real  matrices  or  polynomials,  necessary  and  sufficient 
conditions  have  been  obtained  by  Bialas  [Bl]  and,  in  subsequent  independent  work,  by  Fu 
and  Barmish  [FBI].  Very  recently,  the  case  of  one-  and  two-parameter  families  of  matrices 
was  resolved  for  several  specific  domains  of  the  complex  plane  [GT],  [VT],  [STA]. 

In  this  paper,  we  introduce  and  apply  to  the  generalized  stability  of  parametrized 
families  problem  the  concepts  of  “guardian  map”  and  “semiguardian  map”  for  sets  of 
matrices  or  polynomials.  These  notions  allow  one  to  replace  the  question  at  hand  with 
that  of  whether  or  not  the  guardian  map  is  nonzero  for  all  members  of  the  family  (this 
key  fact  is  proved  in  Section  3).  These  concepts  are  closely  related  to  work  of  Gutman 
[G2]  and  Gutman  and  Jury  [GJ]  on  root  clustering  in  domains  of  the  complex  plane.  The 
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notion  of  “critical  constraints,”  used  in  examples  in  [G2],  is  similar  in  spirit  to  the  idea  of 
guardian  maps. 

Necessary  and  sufficient  conditions  are  given  for  stability  of  one-parameter  families 
relative  to  domains  for  which  the  set  of  generalized  stable  matrices  is  endowed  with  either  a 
guardian  or  semiguardian  polynomial  map.  In  particular,  the  technique  yields  conditions 
for  Hurwitz  and  Schur  (i.e.,  discrete-time)  stability  of  the  convex  hull  of  two  matrices 
or  polynomials.  Combined  with  the  Edge  Theorem,  it  solves  the  problem  of  generalized 
stability  of  arbitrary  poly  topes  of  polynomials. 

For  the  two-parameter  case,  we  consider  stability  of  families  of  matrices  relative  to 
domains  to  which  we  can  associate  a  polynomial  guardian  map.  The  first  step  replaces  the 
two-parameter  problem  by  a  one-parameter  stability  problem  relative  to  a  new  domain. 
The  second  step  employs  a  polynomial  semiguardian  map  associated  with  the  new  domain 
to  obtain  necessary  and  sufficient  conditions  for  stability. 

The  paper  is  organized  as  follows.  Section  2  establishes  notation  and  provides  requisite 
background  material.  The  concepts  of  guardian  and  semiguardian  maps  are  introduced 
in  Section  3.  Two  basic  results  and  several  examples  are  also  presented  there.  These 
results  are  applied  to  one-  and  two-parameter  families  of  matrices  in  Section  4  and  Sec¬ 
tion  5  ,  respectively.  In  Section  6,  a  systematic  procedure  for  constructing  guardian  and 
semiguardian  maps  is  presented  for  problems  of  stability  relative  to  domains  with  polyno¬ 
mial  boundaries.  In  Section  7,  some  applications  of  the  techniques  are  presented.  Finally, 
concluding  remarks  are  given  in  Section  8. 

2.  Preliminaries 

This  section  begins  by  establishing  notation,  and  proceeds  to  a  brief  discussion  of  relevant 
background  material  from  algebra. 

2.1.  Notation 

Arg (s):  Argument  of  the  complex  number  s 

O  O 

d7_  ((C+):  Open  left-half  (right-half)  complex  plane 

T>c:  Complement  of  set  T> 
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T>\  Closure  of  set  T> 
dV:  Boundary  of  set  T> 
int(X>):  Interior  of  set  T> 

Identity  matrix  of  dimension  n  (also  denoted  I  when  n  is  clear  from  the  context) 
A, (A):  Eigenvalue  of  matrix  A 

<r(A):  Eigenvalues  of  matrix  A  (counting  multiplicities) 

Z(p):  Zeros  of  polynomial  p  (counting  multiplicities) 
p':  Derivative  of  polynomial  p 
det(A):  Determinant  of  matrix  A 

0:  Generic  open  subset  of  (L,  symmetric  about  the  real  line 

H:  <F\[0,1] 

0:  <F\  [l,oo) 

Vn:  Set  of  all  real  polynomials  of  degree  at  most  n. 

Sn(V):  Set  of  all  n  x  n  real  matrices  with  spectrum  inside  V  C  (P.  <S„(Q)  is  also 
used  to  denote  the  set  of  all  real  polynomials  of  degree  n  with  zeros 
inside  V.  (Sometimes  denoted  S(T>)  when  n  is  clear  from  the  context.) 

0,  0:  Kronecker  product,  Kronecker  sum 
A  QB:  A  ®  (-B) 

A  ■  B:  Bialternate  product  of  A  and  B  (see  Section  2.2) 

Ai2h  Schlaflian  form  of  order  2  of  matrix  A  (“Upp  er  Schlaflian”;  see  Section  2.2) 

A[ 2] :  Infinitesimal  version  of  A^  (“Lower  Schlaflian”;  see  Section  2.2) 

2.2.  Some  Algebraic  Tools 

Kronecker  product  and  sum 

Given  square  matrices  A  and  B  having  dimension  ni  and  n2,  respectively,  the  Kronecker 
product  (e.g.  [B3] )  of  A  and  B ,  denoted  A  <g)  B,  is  the  square  ni ^-dimensional  matrix 
whose  ijth  ri2  x  ri2  block-entry  is  given  by  aijB.  The  Kronecker  sum  A  ©  B  of  A  and  B 
is  the  niri2-dimensional  matrix  A  ®  J„2  +  Ini  0  B.  Note  that  A  0  A  is  linear  in  A. 

The  eigenvalues  of  A®B  and  A®B  consist  of  the  niri2  products  Xi(A)Xj(B)  and  n\ri2 
sums  A;(A)  +  Xj(B),  respectively,  over  all  ordered  pairs  ( i,j ),  *  =  !,...  ,n i,  ,n 2- 
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In  fact,  this  is  simply  a  special  case  of  the  following  more  general  result.  Let  p  be  a  complex 
polynomial  in  the  variables  xi  and  x2,  given  by 


i+j=N 

p(xl,X2)=  P'JXlX2i  (X) 

i,j=0 

and  consider  the  associated  function  of  two  complex  square  matrices  A  and  B 

i+j—N 

P(A,B):=  V  pijA'  ®  Bj.  (2) 

i,j=0 

Lemma  1.  (Stephanos[S3]).  With  the  notation  above,  the  eigenvalues  of  P(A,B)  consist 
of  the  nin.2  values  p(A;(A),  Xj(B))  over  all  possible  (ordered)  pairs  (i,j),  i  =  l,...,ni,  j  = 
l,...,n2. 


Schlaflian  forms 

The  Schlaflian  forms,2  discussed  next,  have  spectral  properties  akin  to  those  of  the  Kro- 
necker  product  and  sum  with  the  advantage  of  reduced  dimensionality.  Let  x  —  (xj , . . . , 
xn)T  and  p  >  2  be  an  integer.  Denote  by  x^  the  N™ -dimensional  vector  (N™  :=  (n+£  :)) 
formed  by  the  lexicographic  listing  of  all  linearly  independent  terms  of  the  form 


T-Pl  2 
'L  i  •*/2 


n 

xnn>  E^*=P’  Pi  -  °’  *  = 

i=l 


(3) 


For  a  given  n  x  n  matrix  A,  the  associated  (upper)  Schlaflian  matrix  of  order  p  (e.g.,  [A2], 
[B4],  [BZ]),  denoted  is  the  A^-dimensional  square  matrix  defined  by  the  implicit 

relationship 

(Ax)[pl  =  A^xW,  Vx  6  IRn.  (4) 


The  related  form  A[p ]  (the  “lower  Schlaflian  matrix”)  is  defined  as  follows.  Consider 
the  equation  x  =  Ax  for  x  G  IRn  ■  Then  A[p]  is  defined  as  the  coefficient  matrix  in  the 
equation 


dxfpl 

dt 


=  4]xw. 


(5) 


2  These  are  also  referred  to  as  power  transformations  [BZ]. 
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It  can  be  shown  that  A[2]  is  the  directional  derivative  of  along  the  direction  In  [B4] . 
As  such,  A[ 2]  is  linear  in  A. 

The  next  result  is  essentially  the  same  as  results  in  [B4]  and  [BZ]. 

Lemma  2.  The  eigenvalues  of  Ay, j  (resp.  A^)  consist  of  the  N*  sums  (products)  over 
distinct  unordered  multiindices  of  the  form 

Aq  (A)  +  •  •  •  +  \ip  (A)  (resp.  Atl  (A)  x  •  •  •  x  A ip  (A)) .  (6) 


In  contrast,  recall  that  the  eigenvalues  of  the  Kronecker  sum  A  0  A  consist  of  the 
n 2  sums  Ai(A)  +  A j(A)  over  ordered  pairs  (i,j).  In  the  light  of  Lemma  2,  it  is  clear  that 
a(A®  A)  =  <t(A[2]  )  (not  counting  multiplicities).  Hence  the  n^”+1--  x  lower  Schlaflian 

matrix  A[2]  may  be  viewed  as  a  redundancy-free  version  of  the  n 2  x  n2  matrix  A©A,  as  far 
as  the  eigenvalues  are  concerned.  Because  of  this,  A[2]  may  be  used  to  advantage,  instead 
of  A  0  A,  in  the  application  of  some  of  the  results  presented  in  the  sequel.  A  similar 
statement  clearly  holds  for  A^l  vs.  A  ®  A. 

Bialternate  product 

Let  A  and  B  be  n  x  n  matrices.  To  introduce  the  bialternate  product  of  A  and  B,  we 
first  establish  some  notation.  Let  Vn  be  the  ~n(n  —  l)-tuple  consisting  of  pairs  of  integers 
{Pi  q)i  P  =  2, 3, . . . ,  n,  q  =  1, . . .  ,p  —  1,  listed  lexicographically.  That  is, 


=  [(2, 1),  (3, 1),  (3, 2),  (4, 1),  (4, 2),  (4, 3), . ,  (n,  n  -  1)]. 

Denote  by  V"  the  ith  entry  of  Vn.  Denote 


(7) 


f((p,q);(r,s)) 


dps 

bqs 


+  det 


bpr 

dqr 


(8) 


where  the  dependence  of  /  on  A  and  B  is  kept  implicit  for  simplicity.  The  bialternate 
product  of  A  and  B  (e.g.,  [F],  [GJ],  [M]),  denoted  A  -  B,  is  a  |n(n  —  l)-dimensional  square 
matrix  whose  ijtil  entry  is  given  by1 


(A  •  B)y  =  /  (F”;  vp). 


(9) 


1  As  far  as  the  properties  discussed  below  are  concerned,  the  particular  ordering  of  Vn 
is  immaterial.  In  the  literature,  it  is  typically  left  unspecified  ( [F] ,  [G2]  and  [GJ]). 
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Define 


<i(A,A)-.=  Y,'l’riA’'-A’<  (10) 

P>? 

and  denote  the  eigenvalues  of  the  n  x  n  matrix  A  by  Ai, . . . ,  An. 

Lemma  3.  (Stephanos  [S3]).  With  the  notation  above,  the  eigenvalues  of  A)  are  the 
\n(n  —  1)  values 


V*(A i,  -\j)  *  =  2,. . .  ,n  ;  j  =  1, . . .  ,i  -  1. 


(11) 


P,« 


For  example,  if  A  is  3  x  3  then  <r(A  •  A)  =  {AiA2,  Ai  A3,  A2A3}.  In  contrast,  note  that  in 
this  case  cr(A  ®  A)  =  {A2,  Ai  A2,  A1A3,  A2A1;  A2,  A2A3,  A3A1,  A3A2,  A2}.  As  another  example, 
it  is  easily  checked  (e.g.  [G2],  [GJ])  that  if 


Q(A)  =  (A2  -I-  A-  A) 
for  an  n  x  n  matrix  A  with  eigenvalues  A1? . . . ,  An,  then 


(12) 


(A,-A„)2  (A2-A3)2 

2  ’  2 


(A2  -  An)2 

2 


The  Bezoutian 


(13) 


Given  any  polynomial  a(s)  —  ansn  +  •  •  •  +  ais  +  ao,  an  ^  0,  define  the  polynomial  a(s )  := 
s,la(s-1)  =  aosn  +  •  •  •  +  an-i-s  +  an  and  the  matrix 


-  ax 

a2  ■  ■ 

02 

On 

0 

: 

an 

-  On 

0  .  .  . 

0  . 

(14) 


The  Bezoutian  B(a ,  h)  of  two  polynomials  a  and  b  may  then  be  expressed  as  the  n  x  n 
matrix,  n  being  the  largest  of  the  degrees  of  a  and  6,  given  by 

B(a,  b )  :=  S(a)S(b)P  -  S(b)S(a)P,  (15) 
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where  P  is  a  certain  permutation  matrix  [LT]. 

Our  interest  in  the  Bezoutian  stems  from  the  following  result. 

Lemma  4.  The  polynomials  a{s)  and  6(s)  have  no  common  zeros  if  and  only  if  the  asso¬ 
ciated  Bezoutian  B(a,  b )  (or,  equivalently,  the  matrix  S(a)S(b )  —  S(b)S(a ))  is  nonsingular. 

The  same  result  holds  if  instead  of  the  Bezoutian,  the  resultant  (or  Sylvester)  matrix, 
which  has  dimension  2 n,  is  used  [LT]  (see  also  [AS]). 

3.  Guardian  and  Semiguardian  Maps 

In  this  section,  the  concepts  of  guardian  and  semiguardian  maps  are  introduced.  These 
concepts  play  a  key  role  in  subsequent  developments.  Several  examples  are  presented. 

S.l.  Guardian  Maps 

The  definition  of  guardian  maps  given  next  has  been  formulated  with  a  view  toward  prob¬ 
lems  more  general  than  the  generalized  stability  issues  considered  in  this  paper.  (See  e.g. 
Example  3.6  below.)  For  the  purposes  of  this  paper,  the  set  S  of  Definition  1  will  usually 
be  of  the  form  <S(ft),  where 


s(ft)  =  {a  e  mnxn  *(A)  c  ft}  (16) 

for  matrix  stability  problems,  and 

<S(ft)  =  {p  e  Vn  :  Zip)  C  ft)  (17) 

for  polynomial  stability  problems.  Here,  ft  is  an  open  subset  of  the  complex  plane.  Such 
sets  <S(ft)  will  be  referred  to  as  ( generalized )  siabiliiy  sets. 

Definition  1.  Let  X  be  the  set  of  all  n  x  n  square  real  (complex)  matrices,  or  the  set  of 
all  polynomials  of  degree  at  most  n  with  real  (complex)  coefficients,  and  let  S  be  an  open 
subset  of  X.  Let  u  map  X  into  (C.  We  say  that  v  guards  S  if  for  all  x  6  S,  the  equivalence 

x  €  dS  v(x)  =  0  (18) 

holds.  In  this  case,  we  also  say  that  u  is  a  guardian  map  for  S. 
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Proposition  1  below  will  be  used  in  subsequent  sections  to  tackle  the  stability  problem 
for  parametrized  families  of  matrices  or  polynomials  relative  to  domains  of  the  complex 
plane  corresponding  to  guarded  stability  sets.  Let  r  =  (rq, . . .  ,r*)  €  U,  where  U  is  a 
pathwise  connected  subset  of  IRk ,  and  let  x(r)  be  a  matrix  or  polynomial  in  X  which 
depends  continuously  on  the  parameter  vector  r.  Given  an  open  subset  S  of  X ,  we  seek 
basic  conditions  for  x(r )  to  lie  within  S  for  all  values  of  r  in  U. 

Proposition  1.  Let  S  be  guarded  by  the  map  v  and  assume  that  x(r°)  E  S  for  some 
r°  E  U.  Then 

x(r)  E  S  for  all  r  E  U  <=>  v(x(r))  ^  0  for  all  r£[/.  (19) 

Proof.  Suppose  that  r(rJ )  ^  S  for  some  r1  E  U.  By  virtue  of  the  pathwise  connectedness 
of  U,  there  exists  a  curve  { r(t )  :  <  €  [<o5ti]}  within  U,  such  that  r(t o)  =  r°  and  r(ti)  —  r1. 
Now  consider  x(r(t))  as  t  increases  from  to.  Since  x(r°)  E  S,  it  follows  that  there  is  a 
t*  E  (to  ]  such  that 

x(r{t*))  E  dS. 

This  implies  that  there  is  an  r*  E  U  (namely  r*  =  r(t*))  such  that 

x(r*)  E  dS. 

Since  v  guards  S ,  we  conclude  that 


i /(r(r*))  =  0. 


This  proves  sufficiency.  Necessity  follows  from  the  openness  and  guardedness  of  S. 


3.2.  Semiguardian  Maps 


□ 


The  following  generalization  of  the  concept  of  guardian  maps  will  prove  useful  in  the 
development  to  follow. 


Definition  2.  Let  S  and  u  be  as  in  Definition  1.  The  map  v  is  said  to  be  semiguarding 
for  S  if,  for  all  x  E  S,  the  implication 


x  E  dS  =>  u(x)  =  0 


(20) 
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holds.  In  this  case,  we  also  say  that  v  is  a  semiguardian  map  for  S.  An  element  x  £  S  for 
which  u(x)  —  0  is  said  to  be  a  blind  spot  for  (v,S). 

In  the  light  of  Definition  1  and  2,  a  guardian  map  for  a  given  set  S  is  simply  a  semiguardian 
map  for  which  the  corresponding  set  of  blind  spots  is  empty. 

The  next  proposition  is  the  analogue  of  Proposition  1  for  semiguarded  sets  of  matrices 
or  polynomials,  the  main  difference  being  that,  in  this  case,  the  blind  spots  must  be  taken 
into  account.  The  result  will  prove  useful  in  the  study  of  two-parameter  problems  in 
Section  5,  and  is  intended  also  for  sets  for  which  only  a  semiguardian  map  is  known. 

Proposition  2.  Let  S  be  semiguarded  by  v  and  assume  that  x(r° )  £  S  for  some  r°  £  U. 
Then  the  equivalence 

a;(r)  £  S  for  all  r  £  U  x(r)  £  S  for  all  r  £  UCI  (21) 

holds,  where 

UCr  :=  {r  £  U  :  v(x(r))  =  0}.  (22) 

Proof.  Similar  to  that  of  Proposition  1. 

Proposition  2  implies  that  for  the  infinite  family  of  real  matrices  or  polynomials 
(a:(r)  :  r  £  U}  to  be  stable  relative  to  Q,  it  suffices  to  check  that  the  family  {x(r)  :  r  £  UCI} 
is  stable.  In  other  words,  to  establish  that  the  family  {x(r)  :  r  £  17}  is  stable  relative  to 
f2,  one  has  to  ensure  that  u(A(r))  =  0  corresponds  to  the  family  “hitting”  the  blind  spots 
and  not  dS(Q).  In  cases  where  UCI  is  a  finite  set,  Proposition  2  therefore  provides  a  tool 
for  asserting  the  stability  of  the  family  {x(r)  :  r  £  17} . 

Remark  3.1.  The  assumption  x(r°)  £  S  for  some  r°  £  U  appearing  in  Proposition  2  is, 
strictly  speaking,  required  only  in  the  case  UCT  —  0. 

Remark  3.2.  For  some  sets  <S,  the  semiguardian  map  v  factors  as  u  —  u\U2  where  U2(x)  —  0 
implies  x  ^  S  (see  Example  3.9  and  Proposition  7).  In  this  situation,  the  set  UCT  in  the 
Proposition  above  may  be  replaced  by  U*T  :=  {r  £  U  :  rq(x(r))  =  0},  with  the  requirement 
that  the  set  U%r  :=  {r  £  U  :  v2(x(r))  —  0}  be  empty. 

3.3.  Examples 
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We  now  give  some  examples  of  guardian  and  semiguardian  maps  and  the  associated 
guarded  and  semiguarded  sets  of  matrices  or  polynomials.  The  number  of  examples  given 
is  justified  by  the  important  role  each  plays  in  the  development  to  follow.  Before  proceed¬ 
ing  with  the  examples,  we  point  out  that  for  every  guardian  or  semiguardian  map  for  the 
matrix  case  there  are  corresponding  maps  for  the  polynomial  case,  and  vice  versa.  For 
instance,  note  that  a  family  of  polynomials  may  be  viewed  as  characteristic  polynomials 
of  an  associated  family  of  companion  matrices.  Similarly,  stability  of  a  family  of  matrices 
may  be  reformulated  in  terms  of  the  associated  family  of  characteristic  polynomials. 

Example  3.1.  The  map  r:  Ah  det(A)  guards  the  set  of  nonsingular  matrices;  similarly, 
y  :  pH  p(0)  guards  the  set  of  all  polynomials  that  do  not  vanish  at  zero. 

The  next  two  examples  provide  the  simplest  useful  illustrations  of  the  concept  of 
guardian  map:  both  the  set  of  Hurwitz  stable  matrices  (or  polynomials)  and  the  set  of 
Schur  stable  matrices  (or  polynomials)  are  guarded. 

Example  3.2.  The  map  v  :  A  h  det  ( A  ®  A)  guards  the  set  ofnxn  Hurwitz  stable 

O 

matrices  <S((P_).  This  follows  from  the  property  that  the  spectrum  of  A  ®  A  consists  of  all 
pairwise  sums  of  eigenvalues  of  A  (see  Section  2).  Another  such  guardian  map  u  is  given 

O 

by  i /(A)  =  det  (A[2])-  Note  that  each  of  these  maps  guards  5((P+)  as  well.  Similarly,  the 
set  of  Hurwitz  stable  real  polynomials  of  the  form 

p(s)  :=  ansn  +  an-isn~1 - |-ais-fa0  (23) 

is  guarded  by  the  map  i ^  :  pn  det  H(p)  where  H{p)  is  the  Hurwitz  matrix  associated  with 
p  and  is  given  by 


H(P)  = 


On  — 1 

On— 3 

On—  5 

. 

• 

o  - 

On 

On— 2 

On— 4 

. 

0 

0 

On  — 1 

On— 3 

• 

0 

0 

On 

0„-2 

.  . 

• 

0 

0 

. 

•  •  o3 

0 1 

0 

.  0 

. 

. 

02 

Oo  - 

(24) 
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This  follows  from  Orlando’s  formula  [Gl],  [Blj: 


detH(p)  =  a"  Yl  (xi  +  xk), 

'  1  <i<k<n 


(25) 


where  , . . . ,  xn  are  the  zeros  of  p(s ). 

Example  3.3.  The  map  v  :  i  m  det  (A  ®  A  —  I  ®  I)  guards  the  set  of  Schur  stable 
matrices,  i.e.,  of  matrices  with  eigenvalues  in  in  the  open  unit  disk.  Another  guardian  map 
for  this  set  is  given  by  v(A)  =  det(AJ2J  —  I™).  Here  too,  both  maps  guard  the  generalized 
stability  set  corresponding  to  the  outside  of  the  unit  disk.  Similarly,  a  guardian  map  for 
the  set  of  Schur  stable  real  polynomials  of  the  form  (23)  can  be  readily  obtained  from 
results  in  [JP],  [AB]:  Define  the  (n  —  1)  x  (n  —  1)  matrix  D(p)  by 


D{p)  = 


an 

an  —  1 

«n- 2 

.  .  a3 

«2 

-  a0 

0 

an 

an- 1 

.  .  —  do 

-  a\ 

0 

—a  o 

~a\ 

•  •  dfi  — 4 

an— i 

dn_3 

-  —« o 

-02 

•  •  dn_3 

an  _ 

"  dn__2  - 

Then  we  have  the  following  expression  for  det  D(p): 


det D(p)  =  n  (1  -  Xix k). 

k= 1 
«<fc 


(26) 


(27) 


Clearly,  det D(jp)  vanishes  whenever  p  has  a  pair  of  conjuguate  eigenvalues  on  the  unit 
circle.  Note  however  that  detD(p)  does  not  necessarily  vanish  if  p  only  has  1  (or  —  1) 
as  eigenvalue.  Taking  this  into  account,  we  conclude  that  the  set  of  Schur  stable  real 
polynomials  is  guarded  by  u  :  pn  p(l)p(  — l)detD(p). 

In  the  next  two  examples  we  exhibit  guardian  maps  for  sets  of  matrices  having  all 
their  eigenvalues  in  the  domains  depicted  in  Fig.  1  and  Fig.  2. 

Example  S.\.  Given  f)  >  0,  the  map  v  :  A  i— >  det  [(A  +  i/31)  ©  (A  —  i/3I)\  guards  S( 
where  :=  {s  :  |Sm(s)|  <  /?}.  Indeed,  v(A)  =  0  if  and  only  if  A  has  two  eigenvalues 
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Ai  =  x  +  iyi  and  X2  =  x  +  iy2  such  that  y2  -  yi  =  2/3.  Therefore,  if  A  G  S(Q&)  then 
i/( A )  =  0  if  and  only  if  A  has  some  eigenvalues  on  dCl@ ,  i.e.,  A  G  dS( 

Example  3.5.  Given  #o  G  [+ -^,71"),  v  :  d  det  (e’^°  A  ©  e_,s°A)  guards  <S(O^0)  where 

Q,0O  is  the  sector  given  by  Q$0  :=  {s  :  |Arg(s)|  >  0O}.  Again,  i/(A)  =  0  if  and  only  if  A 

has  eigenvalues  A*  =  rkei6k,k  =  1,2  such  that  riei6lei6°  -  r2ei62e~ie°  =  0,  i.e.,  if  and 

only  if  ri  =  r2  and  d2  —  9\  —  29q  (mod  27r).  Since  j  <  6q  <  7r,  this  says  that  for  all 

A  G  S(Qg0),  v(A)  =  0  if  and  only  if  A  has  at  least  one  pair  of  eigenvalues  on  d£lg0  or  a 

0 

single  eigenvalue  at  0.  Note  that  O,6o  =  (C_  for  60  =  j. 

As  a  final  example  of  guardian  maps,  Example  3.6  below  is  of  interest  when  assessing 
“strict  aperiodicity”  of  linear  systems  wherein  the  eigenvalues  (poles)  are  real,  negative 
and  distinct  [J 2] ,  [J3].  It  is  also  a  building  block  for  a  semiguardian  map  to  be  considered 
below. 

Example  3. 6.  The  map  v  :  i  h  det( A2  ■  I  —  A  ■  A)  guards  the  set  of  all  real  matrices 
having  only  algebraically  simple  eigenvalues;  this  follows  directly  from  (11).  Similarly, 
it  follows  from  Lemma  4  that  the  map  u  :  p  >  det B(p,p')  guards  the  set  of  all  real 
polynomials  with  only  algebraically  simple  zeros.  Using  this,  it  can  be  seen  that  the  map 
v  :  A  det(A2  •  I  —  A  ■  A)  det(A)  (resp.  v  :  p  deti?(p,p')  p(0))  guards  the  set  of 
strictly  aperiodic  matrices  (resp.  polynomials). 

Our  first  examples  of  semiguardian  maps  are  related  to  Examples  3.5  and  3.6.  In 
Section  6,  we  shall  provide  a  systematic  procedure  for  constructing  “polynomic”  guardian 
and  semiguardian  maps  for  stability  sets  corresponding  to  a  large  class  of  domains  in  the 
complex  plane  (boundaries  of  which  are  defined  by  polynomials). 

Example  3.7.  Given  f3  >  0  and  O  =  {s  :  |^m(s)|  >  /?},  the  set  Sn(£l)  is  not  guarded 

by  v  :  A  i— >  det  [(A  +  i/3I)  ©  (A  —  i(3I )]  if  n  >  3.  Indeed,  for  n  >  3,  one  can  easily 
construct  a  matrix  A  G  S(Q)  with  a  pair  of  eigenvalues  Ax  =  x  +  *‘yl5  A2  =  x  +  iy2  such 
that  y2  —  yi  =  2/3.  Clearly,  this  cannot  be  done  if  n  <  3,  and  Sn(£l)  is  guarded  by  v  for 
n  =  2,3.  However,  for  any  n,  u  is  a  semiguardian  map  for  Sn (Q)  and  any  matrix  having 
the  eigenvalues  Ai,  A2  above  is  a  blind  spot  for  (v,  Sn(Q)). 

Example  3.8.  For  60  G  (y,7r),  let  the  interior  of  the  complement  of  the  sector  flg0 
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be  denoted  by  ft  =  {s  :  |  Arg  (s)|  <  $o}-  Again  the  set  <Sn(ft)  is  not  guarded  by 
v  :  A  t— ►  det  (et9°  A  ©  e~td°A)  for  any  n  >  2.  Indeed,  if  A  £  <S(ft)  has  eigenvalues 
el9  and  e~'e  where  0  =  n  —  Oo  <  0o,  then  one  eigenvalue  of  el9°  A  ©  e~,9°  A  is  given  by 
e'e°eie  —  e~te°e~l 6  =  etn  —  e~tn  =  0,  implying  v(A)  =  0  although  A  is  a  stable  matrix. 
Here  as  in  the  previous  example,  v  is  a  semiguardian  map  for  <Sn(ft),  for  any  n.  The  set 
of  blind  spots  includes  any  matrix  with  at  least  one  eigenvalue  in  the  mirror  image  (w.r.t. 
the  imaginary  axis)  of  dft  \  {0}  (dashed  lines  in  Fig.  3). 

Our  final  example  will  be  used  in  the  study  of  two-parameter  families  of  matrices  or 
polynomials. 

Example  3.9.  Let  a ,  /3  be  finite  real  numbers,  with  a  ^  /?,  and  let  ft  :=  (T  \  [a,/3].  Then 
the  set  of  stable  matrices  <S(ft)  is  semiguarded  by  the  map 

u:  A  h-*  det  (A2  •  I  —  A  •  A)  det  ((A  -  aI)(A  -  /?/)) .  (28) 

For  ft  =  d7  \  (— oo,  /?],  <S(ft)  is  semiguarded  by  the  map 

v  :  A  ^  det  (A2  •  I  -  A  •  A)  det(A  -  /?!).  (28) 


This  can  be  seen  by  referring  back  to  Example  3.6  and  noting  that  S((T\IR )  is  semiguarded 
by  v  :  A  det  (A2  •  I  —  A  •  A)  .  In  the  case  of  polynomials,  analogous  semiguardian  maps 
are  given  by 

^  pn  detJ5(p,p')  p(a)p(/3).  (30) 


and 

u  :  p  i  *  det B(p,p)  p(/3), 


(31) 


respectively.  Similar  expressions  based  on  the  polynomial  resultant  [LT]  may  also  be  used 
as  semiguardian  maps.  In  all  these  cases,  the  blind  spots  include  any  matrix  (polynomial) 
having  at  least  one  multiple  eigenvalue  (zero)  off  the  real  axis. 

In  the  sequel,  to  highlight  the  stability  aspect  of  this  work  we  choose  to  focus  on  the 
case  of  stability  sets  <S(ft),  with  the  understanding  that  the  results  hold  in  general.  In 
addition,  for  ease  of  exposition  we  discuss  only  the  matrix  case.  The  development  can  be 
carried  out,  mutatis  mutandis ,  for  the  polynomial  case. 
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4.  One-Parameter  Families 


In  this  section,  we  derive  necessary  and  sufficient  conditions  for  stability  of  a  one-parameter 
family  of  matrices 

A(r )  =  Aq  +  rA\  +  •  •  •  +  rmAm  ,  (32) 

r  E  [0,  l],3  relative  to  a  given  domain  C  (P.  In  (32),  Ak,  k  =  1, . . . ,  m,  are  given  n  x  n 
real  matrices.  In  the  remainder  of  this  section,  the  family  (A(r)  :  r  E  [0, 1]}  is  denoted 
by  A. 

All  guardian  and  semiguardian  maps  considered  thus  fax  are  “polynomic,”  in  the  sense 
of  the  following  definition. 

Definition  3.  A  guardian  map  v  is  said  to  be  polynomic  if  it  is  a  polynomial  function  of 
the  entries  (matrix  case)  or  coefficients  (polynomial  case)  of  its  argument. 

The  case  in  which  a  polynomic  guardian  map  for  5(0)  is  available  is  considered  first. 
We  then  consider  the  case  in  which  only  a  polynomic  semiguardian  map  is  available. 

4-1.  Polynomic  Guardian  Maps 

Let  v  be  a  polynomic  guardian  map  for  5(0).  Then  u(A(r))  is  a  polynomial  in  r.  From 
Proposition  1,  it  follows  that  the  family  A  is  stable  relative  to  0  if  and  only  if  (i)  Aq  is 
stable  relative  to  0  (i.e.,  Ao  E  5(0)),  and  (ii)  the  univariate  polynomial  u(A(r))  has  no 
zeros  in  [0,1].  In  other  words,  if  S{ 0)  is  guarded  by  u  and  A(0)  E  5(0)  then  a  necessary 
and  sufficient  condition  for  the  family  A  to  be  stable  relative  to  0  is  that  the  polynomial 
v(A(r))  be  stable  relative  to  S.  In  such  situations,  it  therefore  suffices  to  merely  check 
that  a  scalar  polynomial  has  no  zeros  in  a  certain  interval;  this  can  be  done  using  a  finite 
algorithm  based  on  Sturm  sequences  (see  e.g.  [Jl] ). 

The  necessary  and  sufficient  condition  given  above  assumes  availability  of  a  guardian 
map  v(A)  explicitly  in  the  form  of  a  polynomial  in  the  entries  of  A.  However,  examples 
considered  in  Section  3,  as  well  as  results  on  a  class  of  domains  considered  in  Section  6, 

3  In  fact,  the  discussion  applies  for  r  constrained  to  lie  in  any  interval  (not  necessarily 
compact). 


15 


show  that  guardian  and  semiguardian  maps  often  occur  in  the  form 

u(A)  =  det  F(A)  (33) 

where  JF  is  a  polynomial  mapping  on  IRnXn .  A  necessary  and  sufficient  condition  analogous 
to  the  one  above,  but  not  requiring  expansion  of  the  determinant  (33),  is  now  formulated. 
With  A(r)  as  in  Eq.  (32),  we  may  write 

Jr(A(r))=Y,r‘Fi(Ao,...,Am)  (34) 

i=0 


Note  that 

Fo(Ao, . . . ,  Am)  =  .F(Ao). 


(35) 


In  the  sequel,  F{  denotes  F{  (Ao  , . . . ,  Am )  for  i  —  0, 1, . . . ,  q. 


Theorem  1.  Let  «S(0)  be  guarded  by  a  map  u  of  the  form  (33),  and  let  Ao  G  S(Cl).  Then 
A(r)  G  <S(S7)  for  all  r  G  [0, 1]  if  and  only  if  M(Ao, . . . ,  Am)  G  «S(0)  where  0  =  (F  \  [1,  oo) 
and 


M  (Ao,  . . . ,  Am) 


'  0  /  .  .  0  ' 

.  .  0  I 

_-M0  .  .  .  —Mq-i . 


(36) 


with 

Mi^F-'Fg-i,  t  =  0,...,9-l 


if  q  >  2,  and 

M(Ao,...,Am)  =  -F~1F1 

if  q  =  1. 

Proof.  From  Proposition  1,  we  have  A(r)  £  S(£l)  for  all  r  G  [0, 1]  if  and  only  if 


(37) 

(38) 


i/(A(r))  ^  0  for  all  r  G  [0, 1]. 


(39) 


Since  Ao  G  S(Q)  and  v  guards  S(Q),  it  follows  that  ^(Ao)  =  det.F(Ao)  7^  0.  Therefore  Fo 
is  invertible.  Thus  (33)  and  (34)  imply 

y(A(r))  =  det  F0  det  (/  +  rM,_  1  +  •  •  •  +  rqMo  ).  (40) 
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This  implies  that  v(A(r ))  is  nonvanishing  for  all  r  G  [0, 1]  if  and  only  if 

XO)  :=  det(  fxqI  +  +  . . .  +  /*Ma  +  M0)  ^  0  (41) 

for  all  /x  G  [1,  oo),  where  [x  :=  L.  Since  x(/x)  is  the  characteristic  polynomial  of  M(Ao,  . . . 

>  Am)  if  q  >  1  (also  of  -Mg_i  if  q  =  1),  we  have  that  A(r )  G  5(0)  for  all  r  G  [0, 1],  if 
and  only  if  M(Ao , . . . ,  Am)  has  no  eigenvalues  in  [1,  oo). 

n 

Stability  of  the  matrix  M(Aq,  . . .  ,Am )  relative  to  0  is  therefore  necessary  and  suffi¬ 
cient  for  stability  of  the  family  A  relative  to  the  guarded  set  5(0). 

4-2.  Polynomic  Semiguardian  Maps 

Let  v  be  a  semiguardian  map  for  5(0)  and  assume  that  A(0)  G  5(0).  From  Proposition  2, 
we  have  that  the  family  A  is  stable  relative  to  0  if  and  only  if  A(r)  G  5(0)  for  all  r  G  Ucr 
where  Ucr  is  the  set  of  zeros  of  the  polynomial  u(A(r))  belonging  to  [0, 1]. 

For  semiguardian  maps  of  the  form  (33),  if  the  matrix  M(A0, . . . ,  Am )  is  well-defined, 
the  condition  M(Ao,...,Am)  G  <S(@)  of  Section  4.1  remains  sufficient  for  stability  of 
the  family  A,  but  is  no  longer  necessary.  The  test  is  inconclusive  if  M(A0, . . . ,  Am) 
has  an  eigenvalue  in  [l,oo).  Suppose  that  u(A0)  =  detJF(A0)  ^  0.  Then  the  matrix 
M(Aq  , . . . ,  Am )  given  by  (36)  (or  (38))  is  well  defined.  Define  the  critical  subset  of  the 
spectrum  of  M(A0, . . . ,  Am)  by 

Scr  :=  <j(M(Ao, •  •  • ,  Am))  D  [1,  oo).  (42) 

If  Scr  is  nonempty,  denote  it  by  {/zj, . . . ,  fxe}.  Since  fx  =  i,  the  set  Ucr  =  {r  G  [0, 1]  :  u(A(r)) 
0}  will  then  be  given  by  {jif1,. . .  Proposition  2  now  yields  the  following. 

Theorem  2.  Let  5(0)  be  semiguarded  by  a  map  v  of  the  form  (33),  and  let  u(A0)  ^  0. 
If  Ecr  =  0,  i.e.,  M(A0, . . . ,  Am)  has  no  eigenvalues  in  [1,  oo),  then  the  family  A  is  stable 
relative  to  0  if  and  only  if  Aq  G  5(0).  If,  however,  Scr  =  {fi i, . . .  ,/if}  ^  0,  then  the  family 
A  is  stable  relative  to  O  if  and  only  if 

Aifx-1)  G5(0),  t  =  (43) 
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Remark  4-1-  A  result  analogous  to  Theorem  2  may  be  obtained  with  v  a  polynomial  map 
rather  than  being  specifically  of  the  form  (33).  The  companion  matrix  associated  with 
the  scalar  polynomial  v{A{r))  then  plays  the  role  of  M(Ao,  ■ . . ,  Am),  and  the  assumption 
v(Aq)  7^  0  is  then  no  longer  relevant. 

The  results  above  may  be  applied  to  the  special  case  of  generalized  stability  of  the 
convex  hull  of  two  matrices  or  polynomials.  Combined  with  Edge  Theorem  [BHL],  these 
results  therefore  provide  a  solution  to  the  generalized  stability  of  a  polytope  of  polynomials. 
The  special  case  of  Hurwitz  and  Schur  stability  of  the  convex  hull  of  two  matrices  or 
polynomials  is  presented  next. 

4-3.  Hurwitz  and  Schur  Stability  of  the  Convex  Hull 
of  Two  Matrices  or  Polynomials 

As  applications  of  Theorem  1,  we  consider  Hurwitz  and  Schur  stability  of  the  convex  hull 
of  two  real  matrices.  For  the  former  problem,  a  known  result  is  obtained  [Bl],  [FBI]. 

Given  two  n  x  n  real  matrices  Ao  and  Ai,  the  convex  hull  co(Ao,Ai)  of  Ao  and  Ai 
consists  of  the  matrices 

A(r)  =  (1  —  r)Ao  +  rA\ 

=  A0  +  r{Ai  —  Ao), 

for  r  €  [0, 1]. 

Hurwitz  stability 

O 

Let  Q  =  (F_  and  recall  that  <S(S7)  is  guarded  by  v  :  A  y-y  det^A),  where  F(A)  may 
denote  either  A[2]  or  A  ©  A. 

Corollary  1.  Let  Ao  be  Hurwitz  stable.  Then  co(A0,  A\)  is  Hurwitz  stable  if  and  only  if 
F~1(Aq)T(A\)  has  no  eigenvalues  in  (— oo,0].  Here,  F{A)  can  denote  either  A[2]  or  A®  A. 
Proof.  Since  T  is  linear, 

i/(A(r))  =  det  T{Aq  +  r(Ai  -  A0)) 

=  det  (^(A0)  +  r^(A!  -  Ao)). 

Hence  F0  =  F{Aq)  and  F\  =  F(A\  —  A0),  in  the  notation  of  Theorem  1.  Applying 
Theorem  1  with  q  =  1  yields  that  co(Aq,  Ai)  is  Hurwitz  if  and  only  if  M(A0,  A\  —  Ao)  = 
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—J-  1{Aq)J-’(A\  —  Ao)  G  »S(Q),  i.e.,  has  no  eigenvalue  in  the  interval  [l,oo).  Finally, 
— Fq1F{A\  —  Ao)  =  I  —  F~1{Aq)F{A\),  and  the  result  follows. 

n 

The  analogue  of  Corollary  1  for  the  case  of  the  convex  hull  of  two  polynomials  po  and  p\ 
([Bl,  Theorem  1],  [FBI,  Theorem  3.1])  follows  along  the  same  lines,  from  Example  3.2. 

Corollary  2.  Let  po  be  Hurwitz  stable.  Then  co(po,Pi)  is  Hurwitz  stable  if  and  only  if 
the  n  x  n  matrix  H~1(po)H(p\)  has  no  eigenvalues  in  (— oo,0].  Here,  H(p )  denotes  the 
Hurwitz  matrix  associated  with  the  polynomial  p. 

Discrete-time  (Schur)  stability 

Let  Q.  be  the  open  unit  disk.  From  Example  3.3,  S(Cl)  is  guarded  by  v  :  A  i— ►  det^A) 
where  F(A)  may  be  taken  as  either  A®  A  —  I ®  I  or  A ^  —  I t2l .  Although  the  latter  map 
is  preferable  from  a  computational  point  of  view,  the  former  is  used  here  for  notational 
convenience.  Denoting  A\  —  A0  by  A\ ,  we  have  that 

=  f(A0  +  rAi) 

=  (Ao  ®  Ao)  —  I  <g>  I  +  r  A0  ®  Ai  +  Ai  ®  Ao  +  r2Ai  ®  A\ 

=:  F0 +rFx +r2F2.  (44) 

We  now  apply  Theorem  1  with  F0,FX  and  F2  as  in  Eq.  (44)  and  Ax  identified  with  A\. 

Corollary  3.  Let  all  the  eigenvalues  of  Ao  have  magnitude  less  than  1.  Then  the  same 
is  true  for  any  matrix  in  co(Ao,  Ai)  if  and  only  if  Af(Ao,  Ai)  has  no  eigenvalues  in  [1,  oo), 
where 

M(  A0 ,  Ai )  =  _F-ip2  -F-'F,  (45) 

and  Fq,Fi  and  F2  are  as  in  (44). 

A  similar  result  ([AB])  for  the  case  of  Schur  stability  of  the  convex  hull  of  two  polynomials 
co(Po,Pi)  may  be  obtained  using  the  map  given  in  Example  3.3  and  is  given  next.  For  a 
proof,  see  [Si]. 
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Corollary  4 ■  Let  p0  and  p\  be  Schur  stable.  Then  co(po,pi )  is  Schur  stable  if  and  only 
if  the  matrix  D~1(po)D(pi )  has  no  eigenvalues  in  (  — oo,0).  Here,  the  matrix  D(p)  is  the 
one  defined  by  (26). 

5.  Two-Parameter  Families 

In  this  section,  we  consider  stability  of  two-parameter  families  of  real  matrices  relative 
to  a  domain  for  which  S(Cl)  is  endowed  with  a  polynomic  guardian  map  uq,  which 
we  assume  to  be  real  valued.4  The  matrices  we  study,  denoted  A(ri,r2),  are  polynomic 
functions  of  the  parameters  ri,r2  which  are  taken  to  lie  in  [0,1]. 5  Since  both  uq  and 
A(r  1 ,  r2)  are  polynomic  in  their  arguments,  we  may  express  v$i(A(ri ,  r2))  in  the  form  of  a 
bivariate  polynomial,  i/n(ri,r2),  viz. 

*1  j>2 

vn{rx,r2)  :=  vn{A(rx,r2))  =  ^  Y!  u'i d2r'i  r2  ■  (46) 

ii=0  i2=:0 

Assume  that  A(0,0)  E  S(£l).  Proposition  1  then  implies  that  A(r\,r2)  E  S(Ct)  for  all 
ri,  r2  E  [0,1]  precisely  when 

^a(ri,r2)  ^0  for  all  (rlfr2)  E  [0, 1]  x  [0,1].  (47) 

Note  that  ^n(ri,r2)  is  not  identically  zero,  as  i^0,o  =  ^n(A(0,0))  7^  0  since  A(0,0)  6  S(Q) 
and  since  <S(f2)  is  guarded  by  i/q- 

In  view  of  condition  (47),  the  generalized  stability  problem  being  considered  is  re¬ 
ducible  to  a  positivity  question  (e.g.,  [B2])  which  can  be  addressed  using  work  of  [WZ], 
[AS].  This  approach  was  recently  taken  in  [GT]  and  [VT],  in  addressing  the  two-parameter 
generalized  stability  problem  relative  to  several  specific  domains. 

Notwithstanding  the  amenability  of  the  problem  at  hand  to  solution  by  positivity 
methods,  we  now  consider  this  problem  using  semiguardian  maps.  This  calculation  is  given 

4  All  of  the  examples  considered  thus  far  are  of  this  type.  In  the  cases  of  Examples  3.4 

and  3.5,  real  guardian  maps  can  be  obtained  using  results  of  Section  6. 

5  More  general  intervals  may  be  considered. 
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because  it  leads  one  to  consider  the  issue  of  guardedness  of  approximations  of  5(E)  (only 
a  polynomic  semiguardian  map  is  known  for  5(E)),  which  is  anticipated  to  be  important 
in  multiparameter  problems. 

The  bivariate  polynomial  ,  r2)  can  be  rewritten  in  the  form  of  a  univariate  poly¬ 
nomial  in,  say  r2,  viz. 

*/n(ri,r2)  =  cio(r-i)  +  aq(rq)r2  -+ - b  «a-i (n)r2S2-1  +  as2(n)r2S2  (48) 

where  each  coefficient  o^(ri),  i  =  0, ...,s2  is  a  polynomial  in  r\.  Denote  by  pn  the 
polynomial  in  r2  resulting  upon  fixing  rq  in  (48).  The  bivariate  polynomial  v$i(r\  ,  r2)  does 
not  vanish  for  any  ?q ,  r2  €  [0, 1]  if  and  only  if 

Pn  €  5(E)  for  each  r*i  E  [0,1].  (49) 

Consequently,  a  generalized  stability  question  for  a  two-parameter  family  of  matrices  has 
been  reduced  to  a  similar  question  for  a  related  one-parameter  family  of  polynomials 
relative  to  the  specific  domain  E. 

From  Example  3.9,  the  map  v-~  given  by 

* 'e(p)  =  det B(p,p')  p(0)p(l)  (50) 

is  semiguarding  for  5(E).  Letting  V\  (p),  z/2(p)  denote  det2?(p,p')  andp(0)p(l),  respectively, 
we  obtain  by  combining  Remark  3.2  and  Proposition  2  that  pri  E  5(H)  for  all  rq  E  [0, 1]  if 
and  only  if  the  polynomial  po  G  5(E),  U\x  =  0  and 

Pn  G  5(E)  for  all  rq  G  U£T,  (51) 

where 

Uqi  :=  {ri  G  [0, 1]  :  iq(pn)  =  0},  (52) 

U2cx-.=  {r,  €  [0,1]:  u2(pri)  =  0}.  (53) 

From  (50),  it  is  clear  that  i's(Pn)  is  a  polynomial  in  the  parameter  ?q,  which  we 
assume  not  to  be  identically  zero.  Thus  the  set  UCi  ■—  U^T  U  L,rc2r  is  finite.  For  the  case 
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in  which  UCT  is  empty,  we  have,  by  Remark  3.1,  that  pri  £  5(H)  for  all  r i  £  [0, 1]  if  and 
only  if  po  £  5(H).  Hence,  in  the  current  setting,  A  is  stable  relative  to  0  if  and  only  if 
Po  £  5(H).  Suppose,  on  the  other  hand,  that  U =  0  and  UfT  =:  {^i, . . . ,  pe)  where  the 
Pi's  belong  to  [0, 1].  Then  the  requirement  that  p ^  £  5(H),  i  =  1, . . .  is  necessary  and 
sufficient  for  stability  of  the  family  A  relative  to  0. 

The  following  theorem  summarizes  the  foregoing  discussion. 

Theorem  3.  Let  0  be  a  subset  of  the  complex  plane  such  that  5(0)  is  guarded  by  a  given 
real  polynomic  map  vq.  Let  A  {A(rx,r2)  :  (ri,r2)  £  [oq,/?i]  x  [o:2,/?2]}  be  a  nominally 
stable  family  of  real  matrices;  e.g.,  A(ai,a;2)  £  5(0).  Then  the  family  A  is  stable  relative 
to  O  if  and  only  if  Ufx  =  0  and  the  univariate  polynomials  pai ,  and  pri ,  rq  £  UfT,  have  no 
zeros  in  [a2,/32].  Here,  for  each  ri,  pTx  denotes  the  univariate  polynomial  v^Air i, .)), 

Ulx  :=  {fi  ^  [<*uPi]  ■  det B(pri,p'ri)  =  0} 

and 

■=  {n  e  [<*u0i]  :  PnMPnifo)  =  0}  • 

6.  Techniques  for  Constructing  Guardian  and  Semiguardian  Maps 

6.1.  Generating  New  Guarded  and  Semiguarded  Sets  from  Known  Ones 

The  next  proposition  states  properties  which  provide  means  for  the  construction  of  new 
domains  from  existing  ones  with  the  resulting  generalized  stability  sets  being  guarded  or 
semiguarded. 

Proposition  3.  Let  5,  5i  and  52  be  subsets  of  IRnXn. 

(i)  Assume  that  5  is  guarded  (resp.  semiguarded)  by  u.  Then  —5  :=  {—A  :  A  6  5)  is 
guarded  (resp.  semiguarded)  by  i/_  :  A  >—>  u(  —  A).  In  particular,  if  5  =  5(0),  5(— 0) 
is  guarded  by  i/_,  where  —0  :=  {— s  :  s  £  0}. 

(ii)  Let  5  be  guarded  (resp.  semiguarded)  by  u  and  let  a  £  IR.  Then  5  +  al  := 
{A  +  a  I  :  A  £  5}  is  guarded  (resp.  semiguarded)  by  :  A  i— >  u(A  —  a  I). 
In  particular,  if  5  =  5(0),  5(0*“))  is  guarded  (resp.  semiguarded)  by  Aa\  where 
0(a)  :=  {s  +  a  :  s  £  0}. 
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(iii)  Let  S  be  guarded  (resp.  semiguarded)  by  u  and  let  p  E  IR  \  {0}.  Then  pS  := 
{pA  :  A  €  <9}  is  guarded  (resp.  semiguarded)  by  vp  :  A  t-*  u(~).  In  particular,  if 
S  =  S(Q),  S(pQ)  is  guarded  (resp.  semiguarded)  by  vp. 

(iv)  Let  S\  and  <92  be  guarded  (resp.  semiguarded)  by  v\  and  i/2,  respectively.  Then 
Si  f~l  S2  is  guarded  (resp.  semiguarded)  by  v  :  A  h  i/1(j4)i/2(A).  In  particular,  if 
S\  =  <9(Qi)  and  S2  =  <9(02)5  then  <9(Qi  fl  02)  is  guarded  (resp.  semiguarded)  by  u. 
The  analogous  statement  holds  for  polynomials. 

(v)  Let  S  be  guarded  (resp.  semiguarded)  by  v.  Then  any  connected  component  of  S  is 
guarded  (resp.  semiguarded)  by  v. 

Since  any  convex  domain  (symmetric  w.r.t.  to  the  real  axis,  symmetric  for  short)  with 
polygonal  boundary  may  be  generated  from  the  two  basic  domains  and  0#,  using  the 
basic  operations  in  Proposition  3,  we  have  the  following  subsidiary  result. 

Proposition  4.  Let  Q  be  any  (symmetric)  convex  domain  with  polygonal  boundary. 
Then  <S(fi)  is  guarded  by  a  polynomic  map.  Moreover,  Proposition  3  can  be  used  to 
construct  a  guardian  map. 

6.2.  Domains  with  Polynomial  Boundary 

In  this  section,  we  construct  guardian  and  semiguardian  maps  for  generalized  stability  sets 
corresponding  to  a  whole  class  of  domains  of  the  complex  plane.  Specifically,  we  consider 
domains  whose  boundaries  are  given  by  a  polynomial  equations  p(x,  y)  =  0  where  x  and  y 
denote  real  and  imaginary  parts,  respectively. 

Denote 

D  =  {s  =  x  4-  iy  :  p(x,  y)  <  0} 

where 

P(x,y)  =  ^2 Pktxky2t  ,  (54) 

k,e 

is  a  real  polynomial.  The  fact  that  we  focus  on  real  matrices  is  accounted  for  by  considering 
polynomials  containing  only  even  powers  of  y.  Thus  only  domains  symmetric  w.r.t.  the 
real  axis  are  considered. 
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Associate  with  p  the  real  valued  polynomial 


=  p 


(\  +  \  A  —  A\ 
V  2  ’  2i  ) 


=  1/ 

k,t 


k+ U 

(A  +  A)l(A  -  A)2'. 


(55) 


Rewrite  (55)  as 

q(\,\)  =  J2qkAkXt  (56) 

k,e 

where  the  coefficients  q^t  are  real. 

With  this  notation,  Q  and  dil  have  the  alternative  expressions 


0  =  {A  £  d7 : 

q( A,  A)  <  0} 

(57) 

dn  =  {A  £  <£  : 

?(A,A)  =  0} 

(58) 

Consider  the  mapping  T  :  IRnxn  — >  IRn 

Xn  given  by 

F(A)  :=  ^  qktAk  (g>  A1.  (59) 

k,t 


Lemma  1  implies  that  with  a(A )  =  {A1? . . . ,  A„}, 


<r(^(A))  =  {g(A,-,Aj)  :  i,j  =  (60) 

Now  suppose  that  A  £  dS(Cl).  Then  some  eigenvalue  of  A  satisfies  A;  £  dQ,  i.e.,  q(\i ,  A;)  = 
0.  It  then  follows  from  (60)  that  J~{A)  is  singular  (det  J~{A)  =  0).  We  obtain  the  following 
propositions. 

Proposition  5.  Assume  that  v  is  not  identically  zero.  Then  the  map 


v  :  A  t-»  det  <lktAk  ®  A e 
k,e 


is  semiguarding  for  <5(0). 


(61) 
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Proposition  6.  The  map  (61)  guards  5(0)  if  and  only  if  q  satisfies  the  condition7 


#(A,A)<0  and  <?(//,//)<  0  =>  q(\,fi)^0. 


(Condition  C ) 


Maps  such  as  (61)  involve  determinants  of  matrices  the  size  of  which  increases  rapidly  as  n 
does.  Alternative  formulas,  based  on  the  bialternate  product,  exist  which  involve  matrices 
of  dimension  — ,  which  is  approximately  half  that  of  J-{A)  for  large  n. 


Consider  the  mapping  A4  from  IRn  n  to  IR 


(  n(n-l)  x  n(n-l) 


2  given  by 


M(A)  =  ]TqksAk-A*. 
M 


(62) 


Lemma  3  implies  that 


cr(M(A))  =  {X><(A*A*  +  AfAj)/2  :  i  =  1, . 
^  k,e 


, n  -  1;  j  =  i  +  1,...  ,n 


(63) 


Note  that 

<r(Ad(A))  =  {<?(Ai,Aj>  :  i  =  1, . . . ,  n  -  1;  j  =  i  +  1, . . . ,  n}.  (64) 

(Compare  with  Eq.  (60).)  This  follows  from  (63)  and  the  fact  (implied  by  (55))  that 

q( A,  ft)  =  q(n,  A)  =  ^  (g(A,  //)  +  q(n,  A)) .  (65) 

Proposition  7.  Suppose  that  detAt(A)  is  not  identically  zero. 

(a)  If  <90  fl  iR  =  0,  then  the  map 

v  :  A  i— »  detA4(A)  (66) 

is  semiguarding  for  5(0). 

7  “O-transformability”  [G J] . 
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(b)  Let  dH  {MR  =  [a*,  /3,-j  with  oq  <  P\  <•■■<  ai  <  Pi  <■••<  <  P$.  Denote 

ha)  =  nf=.(^  —  a.iI)(A  —  flil)  where,  by  convention,  the  factor  ( A  —  oq/ )  (resp. 
( A  —  P^I))  is  omitted  when  aq  (resp.  (3^)  is  — oo  (resp.  +oo).  Then  the  mapping 

v  :  Am  detAf (A)det'P(A)  (67) 

is  semiguarding  for  S(Q). 

Proof.  Let  A  E  IRnXn  fl  dS(Q): 

(a)  If  A  has  an  eigenvalue  A  £  dtt  then  A,  also  an  eigenvalue  of  A ,  is  distinct  from  A. 
Therefore  ^(A,  A)  =  0  £  cr(M(A))  by  virtue  of  (62).  Hence  u(A)  =  0. 

(b)  Let  A  £  dQ  be  an  eigenvalue  of  A.  That  is,  q( A,  A)  =  0.  If  A  ^  IR,  then  A  is  also  an 
eigenvalue  of  A,  distinct  from  A.  Consequently,  q( A,  A)  =  0  €  cr(M(A))  and  u(A)  =  0. 
If  A  £  IRndQ,,  then  A  £  [a,-,  /?*•]  for  some  i  £  {1, . . . ,  £}.  Since  by  assumption  A  E  d£l, 
then  it  is  the  limit  of  a  sequence  of  matrices  {^4fc}  with  each  Ak  E  S(CL).  It  follows 
that  there  is  a  j  E  {l,...n}  such  that  A  =  limfc—oo  Xj(Ak).  If  A  E  (oq,/?i)  (Pi  >  oq 
assumed)  then  there  is  a  positive  K  such  that  for  all  k  >  K 

^j(Ak)  £  (U  \  IR  and  A=  lim  Xj(Ak). 

k—*  oo 

Since  is  a  sequence  of  real  matrices,  we  have  that  A  =  limjt_.00  Xj(Ak)  as  well. 

Consequently,  A  must  be  an  eigenvalue  of  A  of  multiplicity  at  least  2.  By  virtue  of 
(62),  q( A,  A)  =  q( A,  A)  E  cr(M(A)),  i.e.,  v(A)  =  0.  If  a*  ^  p,  and  A  is  either  aq  or 
Pi,  then  A  might  be  a  simple  eigenvalue  of  A ,  in  which  case  q( A,  A)  =  0  is  no  longer 
an  eigenvalue  of  M(A).  The  case  when  a,-  and  Pi  are  finite  and  A  =  oq  =  Pi  is 
handled  similarly.  The  reason  for  the  introduction  of  the  second  factor  det'P(AL)  in 
the  expression  (67)  for  u  should  now  be  clear. 


□ 

Remark  6.1.  Proposition  6  applies  for  the  maps  of  Proposition  7  as  well. 

We  have  exhibited  semiguardian  maps  for  generalized  stability  sets  corresponding  to 
domains  with  polynomial  boundary.  Determining  whether  or  not  these  maps  are  also 
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guardian  maps  requires  further  investigation.  One  needs  to  check  whether  or  not  the 
polynomial  q  satisfies  Condition  C.  Sufficient  conditions  for  Condition  C  were  obtained  by 
Gutman  and  Jury  [GJ]  for  the  cases  in  which  the  degree  of  polynomial  p  (or  q)  is  1, 2, 3  or 
4.  Another  result  in  this  direction  is  given  next. 

Proposition  8.  Suppose  that 

qkk  >  0,  V  k  >  1,  (68) 

qkt  =  0,  K^O.  (69) 


Then  q  satisfies  Condition  C. 


Proof.  Proceeding  by  contradiction,  assume  that  for  some  pair  (A,//),  q( p)  =  0, 
q( A,  A)  <  0  and  q(p,p)  <  0.  Since  the  coefficients  of  q  are  real,  we  also  have  q(X,p)  =  0. 
Set 


w  :=q(X,X)  +  q(p,p). 


Clearly  w  <  0  and 


W  =  q{ A,  A)  +  q(p,  p)  -  (q( A,  p)  +  q( A,  p)) 

=  qke(Xk\e  +  pk  pe  —  (A k  pe  +  A  kpe)). 

k,t 


From  (55)  and  (56),  we  have 

qko  -  qok,  Vfc  >  1.  (70) 

It  now  follows  from  (69)  and  (70)  that 

w  =  ]T  (Xk  +  pk-  (Xk  +  A*)  +  A k+pk-  ( pk  +  pk )) 

k=l 

+  J2 «*»  (A‘A* + /<  v  -  <A  v + a1 p i)). 

fc=i 

The  first  summation  yields  zero.  Under  assumption  (68),  the  second  summation  is  non¬ 
negative,  as  can  be  seen  by  noting  that,  with  a  =  r\e^x  and  (3  =  r2e^2 , 

oca  +  (3/3  —  (a [3  +  a (3)  —  r\  +  r%  —  2r\r2  cos  {8\  +  d2) 

>  ( ri  -  r2)2. 
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This  contradicts  the  fact  that  w  <  0. 


□ 

Example  6.1.  Let  Q,  =  {.s  =  x  +  iy  :  x  -f-  y2  <  0}.  Here  p(x,y)  =  x  +  y2  and 

It  follows  from  Proposition  8  that  q  satisfies  Condition  C,  and  therefore  that  S(Cl)  is 
guarded  by  both 


v  :  A  i — ►  det^A)  and  v  :  A  t— >  det  A^(^4)det(^4). 

Here, 

X(A)  =  t(A®  I+I®A)  +  lA®A-hA2®I+I®A2) 

=  l(AS>A)-l(AQAf,  (71) 

and 

M{A)  =  -(A  ■  I  +  I  ■  A)  +  -A  •  A  -  -{ A 2  ■  I  +  I  ■  A2) 

2  2  4 

=  AI+^A-A-ll-A2.  (72) 

7.  Examples  of  Application 

We  apply  some  of  the  results  obtained  in  this  paper  to  three  examples.  The  first 
two  examples  deal  with  the  stability  of  one-parameter  families  of  matrices  relative  to 
certain  domains  which  arise  in  control  system  design.  For  the  second  of  these  examples, 
we  compute  the  largest  interval  of  parameter  variation  containing  0  for  which  stability  is 
preserved.  The  final  example  concerns  the  Hurwitz  stability  of  a  two-parameter  family 
of  matrices.  Programs  implementing  the  techniques  in  this  paper  and  using  the  symbolic 
language  MACSYMA,8  were  developed  and  applied  to  construct  the  examples  [S2]. 

8  MACSYMA  is  a  registered  trademark  of  Symbolics,  Inc.,  Cambridge,  MA. 
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7.1.  Robust  Stability  with  Adequate  Step  Response. 

To  ensure  adequate  maximum  settling  time  and  minimum  damping,  it  is  often  required 
in  the  design  of  compensators  for  linear  control  systems  that  the  eigenvalues  (poles)  of 
the  closed  loop  system  be  confined  within  a  domain  Cl  :=  {s  :  3Je(s)  <  —  a  ;  |Arg(s)|  > 
9  },  for  given  a  >  0  and  9  6  (f,7r)  (see  Fig.  4)  [V].  In  this  example  we  investigate  the 
stability,  with  respect  to  Cl,  of  the  one-parameter  family  of  matrices  A  given  by 


A(r) 


r  —  3  1  2r  +  1 

r  —1  —1 

1  r  +  1  -3 


r  1 

r  E  ^  2’ 2 


5 


in  the  case  a  =  —1  and  9  —  We  first  provide  a  guardian  map  for  <S(fi).  Letting  fig 
denote  the  set  of  Example  3.5,  and  using  the  notation  of  Proposition  3,  we  observe  that 

o  (  —  <r)  o 

Cl  —  Qq  fl  <S _  .  Since  guardian  maps  for  S(Clo)  and  S(d7_)  are  available,  it  follows  from 

Proposition  3  that  S(Cl)  is  guarded  by  the  map  v  specified  by 


v(A)  —  det(e*0A©e  lB A)  det((A  +  <r/)[2])- 

For  a  3  x  3  matrix, 


2an 

2ai2 

2013 

0 

0 

0 

021 

Q11  +  «22 

023 

Oi2 

013 

0 

«31 

O32 

On  +  O33 

0 

012 

O13 

0 

2o2i 

0 

2a22 

2023 

0 

0 

031 

021 

032 

022  +  033 

«23 

0 

0 

2o3i 

0 

2032 

2033 

(73) 


(74) 


Let  us  for  convenience  write  v(A)  =  ui(A)u2(A) 


where  the  z';’s  are  the  respective  factors 


appearing  in  (73).  We  obtain 


vi(A(r))  =  iV2(-32r9  -  16r8  -  344r7  +  336r6  +  663 6r5  +  13604r4  +  15408r3 

-  48008r2  +  38298r  -  22536),  (75) 

u2(A(r))  =  -32 r6  -  192r5  +  288r3  +  64r2  -  640r  +  312.  (76) 


It  is  easily  checked  that  A(0)  G.  S(Cl ),  and,  using  Sturm  sequences  for  example,  that  v\  and 
v2  have  no  zeros  in  [—f ,  f  ]•  We  may  conclude  by  virtue  of  Proposition  1  that  the  family 
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A  is  stable  relative  to  12.  Parenthetically,  the  real  zeros  of  of  v\  and  u2  are  given  by 


{-5.7588,0.5919}  and  {-2.9859,0.95126,3.737}, 


respectively.  Thus,  an  eigenvalue  of  A(r)  leaves  12,  for  the  “first  time,”  at  r  =  0.5919. 

7.5.  Largest  Interval  for  Generalized  Stability 
In  this  example,  we  investigate  the  stability  of  the  one-parameter  family  of  matrices 

r2  —  1  r  +  1 
r2  —  2r  —  1  —1 


relative  to  the  domain  (see  Fig.  5)  12  =  {s  :  pi  <  |s|  <  p2  ;  |Arg(s)|  >  #},  with 


Pi 


V2 

2  ’ 


P2  =  2. 


Such  domains  arise  in  flight  controller  design  [Al,  p.  394].  We  seek  to  compute  the  largest 
open  interval  (r,r),  containing  0,  such  that  A(r)  G  «S(0)  for  all  r  G  (r,r)  (note  that  the 
eigenvalues  of  A(0)  are  —1  ±  t,  hence  A(0)  G  <S(12)).  To  apply  Proposition  3,  we  write 


12  =  12#  fl  B(p\)  fl  int(Bc(p2)) 

where  B(p)  denotes  the  open  disk  of  radius  p.  Thus  «S’(122 )  is  guarded  by  the  map  v2  given 
by 

v2{A)  =  det (eidA  0  e~i$A )  det(^2I  -  p\l[2])  det(A^  -  p2I[2]) 

=:  vi(A)u2(A)v3(A). 


Here,  A^  is  given  by 


a2j  2anai2  a22 

ana2l  0-110,22  +  a12a21  O12O22 

2  o  2 

O21  2021^22  °22 

For  the  example  at  hand  we  obtain 

i/i(4(r))  =  — 3r7  -  3r6  +  21r5  +  24r4  -  36r3  -  51r2  +  12, 

v2(A{r))  =  l  (8r9  -  76r7  -  52r6  +  244r5  +  318r4  -  170r3  -  492r2  -  282 r  -  5l) 
8 

va(A(r))  =  -r9  +  13 r7  +  10r6  -  55r5  -  52r4  +  79r3  +  58r2  -  12r  -  40. 


(77) 
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The  real  zeros  of  the  polynomials  above  are 


r  —  —2,  r  ~  —1.246,  r  =  —1,  r  ~  0.445,  r  ~  1.801,  r  =  2; 

r  «  -1.384,  r  «  -1.112,  r  «  -0.5578,  r  «  -0.4332,  r  »  1.919,  r  «  1.942,  r  «  2.252; 

r  =  —2,  ,  r  =  1,  r  =  2,  r  «  3.152; 

respectively.  Clearly,  the  matrix  A(r)  is  not  stable  relative  to  Cl  for  all  values  of  r  6  77?.. 
It  is  also  seen  furthermore  that  the  maximal  interval  of  stability  is  [—0.4332,  0.445]. 

7.S.  Eurwitz  Stability  of  a  Two -Parameter  Family 

O 

In  this  example,  Cl  =  (T_  and 


A(n,r2) 


— 3  —  r2  +  3t*2  —  l  +  r2+4rir2 

-l  +  2ri  -2  +  3  rx+r2-r|  ’ 


ru  r2  £  [0, 1].  Note  that  A(0,0)  is  Hurwitz  stable.  A  guardian  map  for  is  given  by 


vsi(A(ri,r2))  =  deti4[2](ri,r2), 


where  A^  is  given  by 

2an  2a2i  0 

a21  on  +  a  22  a12 

0  2ai2  2a22 


(78) 


We  obtain 

i 'n(4(r1}r2))  =  -(100  -  200ri  +  84r?)  +  (20rx  +  148r?  -  9 6r*)r2  +  (120  -  284rx  +  100 rj)r22 
-  (80  -  40rx  +  64 r?)rf  +  (28  +  36 rx)r\  +  32 r\  -  24 r® 

=:  Pn(r2). 

In  this  case,  we  find  that 

Ult  «  {0.5689,  0.5725,  0.8207}  and  U c2r  «  {0.5663,  0.7142,  1}  ^  0. 


It  follows  from  Theorem  3  that  the  family  of  matrices  being  considered  is  not  Hurwitz 
stable.  For  rx  =  0.5663,  pri  has  (r2  =)  1.0  as  zero,  implying  that  A(0.5663, 1.0)  is  Hurwitz 
unstable  (its  eigenvalues  are  {  —  1.3, 0.0}). 
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8.  Concluding  Remarks 


This  paper  has  introduced  a  new  framework  for  the  study  of  generalized  stability  of 
parametrized  families  of  matrices  and  polynomials.  Stability  relative  to  a  wide  variety 
of  domains  in  the  complex  plane  can  be  addressed  within  this  framework. 

Proposition  1  and  2  give  necessary  and  sufficient  conditions  for  generalized  stability 
of  multiparameter  problems.  Where  Proposition  1  applies  and  the  family  of  interest  is 
polynomic,  the  problem  is  reducible  to  a  positivity  test. 

Although  the  presentation  has  emphasized  the  case  of  generalized  stability  sets  S(Q), 
the  results  apply  as  well  to  other  open  sets  of  matrices  or  polynomials.  A  case  in  point  is 
the  set  of  strictly  aperiodic  matrices  (see  Example  3.6). 

Finally,  the  results  presented  in  this  paper  can  be  extended  in  a  straightforward  man¬ 
ner  to  the  case  of  matrices  with  complex  entries  or  polynomials  with  complex  coefficients. 
In  particular,  guardian  and  semiguardian  maps  can  be  readily  constructed  for  sets  of  com¬ 
plex  matrices  corresponding  to  domains  of  the  complex  plane  with  polynomial  boundaries. 
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